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Abstract 

A new web-server called The Molecule Calculator (MolCalc) is pre- 
sented. The entry page is a molecular editor (JSmol) for interactive 
molecule building. The resulting structure can then be used to esti- 
mate molecular properties such as heats of formation and other ther- 
modynamic properties, vibrational frequencies and vibrational modes, 
and molecular orbitals and orbital energies. These properties are com- 
puted using the GAMESS program at either the RHF/STO-3G (or- 
bitals and orbital energies) or PM3 level of theory (all other proper- 
ties) in a matter of seconds or minutes depending on the size of the 
molecule. The results, though approximate, can help students develop 
a "chemical intuition" about how molecular structure affects molecu- 
lar properties, without performing the underlying calculations by hand 
- a near impossible task for all but the simplest chemical systems. 
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Introduction 

The Molecule Calculator (MolCalc) is a web interface that allows students 
to build small molecules and estimate molecular properties such as molecu- 
lar structure, thermodynamic properties, vibrational frequencies, vibrational 
modes, molecular orbitals and orbital energies in a matter of seconds or min- 
utes - depending on the size of the molecule. A video demonstrating these 
features can be be found on the website. ® 

MolCalc is designed for teaching and assignments in which the students build 
their own molecules and estimates their own molecular properties, as opposed 
to reading from tables in textbooks. MolCalc is therefore designed to run fast 
and the estimated molecular properties will not match experimental values 
exactly, and in some cases be quite different. The idea is to have students de- 
velop a "chemical intuition" about how molecular structure affects molecular 
properties, without performing the underlying calculations by hand (which 
would be near impossible for all but the simplest chemical systems). 

Just like a pocket calculator or a symbolic math program (such as Math- 
ematica or MAPLE), MolCalc allows one to assign "higher level" chemical 
problems that are not practically possible to solve otherwise. For example, 
one might now ask students to compute the effect of a substituent on a par- 
ticular vibration and then rationalize the effect using molecular orbitals. Or 
one might ask more open ended questions such as "build a molecule with an 
unusually long C-C single bond" . ® 

How MolCalc works 

MolCalc is split into two steps, building the molecule and calculating prop- 
erties. The editor is centered around the JSmol® plugin which is a new 
JavaScript version of the popular Jmol® visualization program that, unlike 
Jmol, works on a broad range of devices without the use af Java, such as the 
iPhone and iPad. 
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Different starting structures are provided, as well as the possibility to search 
for other structures using common names or SMILES. The editor interface is 
based on a point-and-click approach with a clear intuitive workflow, as seen 
in Figure [TJ JSmol is used for all molecule modifications as well as structure 
minimization using the MMFF force field®. 

In the Molecule Calculator page the structure is re-optimized at the PM3 level 
of theory® for a maximum of 50 steps, using the GAMESS program.^ 
After the minimisation the student is presented with the PM3 minimized 
structure and the different calculation types. The PM3 structure can then 
used to compute the heat of formation and other thermodynamic proper- 
ties, vibrational frequencies and normal modes at the PM3 level of theory, or 
the molecular orbitals and orbital energies using the RHF/STO-3G level of 
theory. These calculations are also performed with the GAMESS program. 
OpenBabel® is used to manage input files and coordinate files. 

In order to ensure a fast turnaround MolCalc 1.1 allows calculations on 
(closed shell) molecules with only doubly occupied molecular orbitals and 
with less than 11 non-hydrogen atoms. Only elements up to Ar are al- 
lowed. PM3 is chosen for its speed and the fact that it provides heats of 
formation. RHF/STO-3G is chosen for the orbital calculations because it 
provides the virtual orbitals, particularly the lowest unoccupied molecular 
orbital (LUMO) expected from the linear combination of atomic orbitals 
(LCAO) model. This is not always the case for larger basis sets. 

MolCalc is distributed through Github under the GPL license. ® The Github 
page also accepts bug reports and feature requests. The core of MolCalc is 
written in PHP5, j Query, HTML5 and Python2.7, but is very modular and 
can be installed on most webservers, which makes it quite easy to add new 
capabilities to MolCalc. A copy of the GAMESS code must be obtained 
separately separately. ® 
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Figure 1: The Molecule Calculator (MolCalc) editor page. 



Acknowledgements 



The development of MolCalc is supported by the University of Copenhagen 
through the Education at its Best initiative (Den gode uddannelse). The 
Molecule Editor was inspired, in part, by the Virtual Molecular Modeling 
Kit® developed by Othis Rothenberger and Thomas Newton, and based on 
a prototype developed by Toke Fritzemeier. 



Author Information 

Corresponding Author 

*E-mail: jhjensen@chem.ku.dk 
Notes 

The authors declare no competing financial interest. 

References 



[1] http://dgu.ki.ku.dk/molcalc (accessed March 2013) 



[2] http : / / dgu .ki.ku.dk/ molcalc / about | (accessed March 2013) 



[3] http://github.com/jensengroup/molcalc (accessed March 2013) 



[4] http://molecularmodelingbasics.blogspot.dk/2011/12/bond-lengths- 



avogadro-google -docs-and.html| 



[5] http://jmol.sourceforge.net (accessed March 2013) 



[6] |http: / / che mapps.stolaf.edu / j mol / j smol /test 2 . htm (accessed March 
2013) 

[7] Halgren, T. A., Merck molecular force field. I. Basis, form, scope, pa- 
rameterization, and performance of MMFF94. J. Comp. Chem., 1996, 
17: 490-519. 



[8] http://www.msg.ameslab.gov/gamess/download.html (accessed March 



2013) 

Stewart J. J. P., Optimization of Parameters for Semiempirical Methods 
I. Method J. Comp. Chem. 1989 V. 10 2:209-220 



5 



[10] O'Boyle N. M., Banck M., James C. A., Morley C, Vandermeersch 
T. and Hutchison G. R., Open Babel: An open chemical toolbox J. 
Cheminformatics 2011, 3:33 

[11] Schmidt M.W., Baldridge K.K., Boatz J.A., Elbert S.T., Gordon M.S., 
Jensen J.H., Koseki S., Matsunaga N., Nguyen K.A., Su S., Windus 
T.L., Dupuis M., Montgomery J. A., General Atomic and Molecular Elec- 
tronic Structure System, J. Comput. Chem., 1993, 14, 1347-1363. 

[12] |http://mol ecularmodelingbasics.blogspot.dk/2012/04/virtual- 
molecular-modeling-kit.html (accessed March 2013) 



6 



